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ABSTRACT 



o 
u 

O We present here photometric redshift confirmation of the presence of large scale 

^ structure around the z = 1.82 QSO RX J0941, which shows an overdensity of submm 

^ sources. Radio imaging confirms the presence of the submm sources and pinpoints 

I ^ I their likely optical-NIR counterparts. Four of the five submm sources present in this 

field (including the QSO) have counterparts with redshifts compatible with z = 1.82. 
T-H We show that our photometric redshifts are robust against the use of different spectral 

^ templates. We have measured the galaxy stellar mass of the submm galaxies from their 

rest-frame iC-band luminosity obtaining log(M*/M0) ^ 11.5±0.2, slightly larger than 
the Schechter mass of present day galaxies, and hence indicating that most of the stel- 
lar mass is already formed. We present optical-to-radio spectral energy distributions 
(SEDs) of the five SCUBA sources. The emission of RX J0941 is dominated by re- 
7—i processed AGN emission in the observed MIR range, while the starburst contribution 

completely dominates in the submm range. The SEDs of the other three counterparts 
are compatible with a dominant starburst contribution above ~ 24 /urn, with star 
formation rates SFR^ 2OOOM0/yr, central dust masses log(Mdust/-^0) ^ 9 ± 0.5 
and hence central gas masses log(Mgas/M0) ~ 10.7. There is very little room for an 
AGN contribution. From X-ray upper limits and the observed 24 /xm flux, we derive 
a maximum 2-lOkeV X-ray luminosity of 10''^ erg/s for any putative AGN, even if 
they are heavily obscured. This in turn points to relatively small black holes with 
log(M,/M0) < 8 and hence stellar-to-black hole mass ratios about one order of mag- 
nitude higher than those observed in the present Universe: most of their central black 
hole masses are still to be accreted. Local stellar-to-black hole masses ratios can be 
reached if ^ 1.3% of the available nuclear gas mass is accreted. 

Key words: galaxies: formation - galaxies: starburst - galaxies: evolution - galaxies: 
high redshift - submilimetre 
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1 INTRODUCTION 

It is currently commonly accepted that structure in the 
Universe is formed through hierarchical processes, in which 
smaller structures collapse to form galaxies, which then 
group together or coalesce to form larger galaxies and galaxy 
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groups and clusters. Within this general background, "anti- 
hierarchical" behaviour has been observed, in the sense of 
the most dense large-galaxy-size structures collapsing ear- 



lier and evolving faster than smaller galaxies (e.g. Renzini 
2006 and references therein). These processes are thought 



to be accompanied by channeling of material to the central 
regions, where it form stars and feeds a growing black hole 
(BH). The feedback of the latter regulates galaxy formation 
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in an "evolutionary sequence" (e.g. Silk & Rees|1998 


Fabian 


19991 IGranato et al.||2004| |Page et al.||2004 King||2005 Di 


Matteo, Springel, & Hernquist"2005| 


Stevens et al. 


2005 


1. 



According to these models galaxies grow through star for- 
mation and at first host "small" BH possibly in very dense 
and obscured nuclear environments. As time progresses the 
combined radiative and mechanical output from star forma- 
tion and the BIf accretion grows too, until the AGN reaches 
QSO luminosities and literally blows away the circumnuclear 
material. This effectively terminates star formation and the 
SMBH shines briefly (in cosmic terms) as an unobscured 
QSO. Once it has accreted the remaining material, the QSO 
switches off, leaving a passively evolving galaxy with a "dor- 
mant" SMBff, such as those observed in present-day massive 
galaxies ( [Marconi fc Ifunt||2003 1. 



The Far-infrared (FIR)-submillimetric (submm) spec- 
tral region is ideally suited for searching for star formation, 
since the black-body-like (greybody) emission by cold dust 
associated with star formation peaks at those wavelengths. 
At the same time it suffers very modest obscuration from the 
circumnuclear material. This is proved by the large numbers 
of strongly star forming submm galaxies (SMGs) detected by 
SCUBA, AzTEC, SMA and, lately, Herschel, among other 
facilities. 

There are many multi-wavelength studies of SMGs, 



both from surveys in "blank fields" (e.g. |Smail, Ivison, & 


Blain 


19971 Hughes et al.ll998J Bales et al.|1999| Scott et al. 


2002 


Borys et al.|2005, Alexander et al.|2005| Coppin et al. 


2006 


Laird et al.||2010) and targeting overdensities of such 



sources around high rcdshift objects (e.g.,|Kurk et al.|2000 



Pentericci et al. 2000; Ivison ct al. 2000; Stevens ct al. 2003 

2004,.Greve et al..,2004 



^maiT^ al.||2003j ,De Breuck et al 

Stevens et al.|2004[|Venemans et al.|2007|rPriddey, Ivison, fc| 



Isaak||20'08 Chapman et al.||2009 K These studies probe the 
properties of the star formation and of the AGN, and the as- 
sociation between both phenomena, although there are not 
many cases in which this association has been proved. SMGs 
have substantial stellar masses in place by z ~ 2. A relatively 
high fraction of SMGs host AGN (e.g. [Alexander et al.|2005 



but see Laird et al.|2010 \ , which seem to have BH-to-stellar 



mass ratios smaller than local galaxies ( jBorys et al.[|2005| 



Alexander et al. 20081. This implies that BH growth lags 
galaxy growth, since the host galaxies are already mature, 
while the BHs still require substantial growth (e.g., by about 
a factor of 6), which however can be accommodated within 
the limits imposed by the lifetime of the submm-bright phase 
( Alexander et al.|2008 l, if accretion occurs close to the Ed- 
dington limit ( [Eddington 1913 Rees 19841. The observed 
population of z ~ 2 SMGs would be sufficient to account for 
the formation of the population of bright ellipticals seen in 



the present universe ( Swinbank et al. 2006 1 



We have found ( [Page et al.||2004| [Stevens et al.||2005| ) 
a sample of X-ray-obscured QSOs at 2; ~ 2 (when most 
of the star formation and BH growth is occuring in the 
Universe) with strong submm emission, much higher than 
X-ray-unobscured QSOs at similar redshifs and luminosi- 
ties, which however represent 85-90% of the X-ray QSOs at 
that epoch. This submm emission is too high to come from 
the AGN ( [Page et al.[|2001[ [Stevens et al.|[2005[ ). The fields 
around these objects show strong overdensities of SMGs 
( [Stevens et al.|2004[[20Tol). The UV and X-ray spectra of the 
central QSOs (Page et al.|2010 ) show evidence for strong ion- 



ized winds which produce the X-ray obscuration. Piecing all 
these clues together we infer that the host galaxies of these 
QSOs are undergoing strong SFR, while the central SMBH 
are also growing through accretion. The ionized winds are 
strong enough to quench the star formation, so the QSOs 
must be just emerging from a strongly obscured accretion 
state to become "normal" unobscured QSOs with passively 
evolving galaxies, in an evolutionary stage which might last 
about 10-15% of the QSO lifetime. They appear to be in the 
centres of high density peaks of the Universe. 

In this paper we endeavour to prove the physical as- 
sociation of one of those overdensities to its central QSO 
(RX J094144.51-f-385434.8, henceforth RX J0941). We also 
investigate the nature and evolutionary stage of those SMGs 
through their rest frame UV-FIR SEDs. 

The outline of this paper is as follows: in Section [2] we 
summarize the data used in this paper and the reduction 
process, in Section [3] we describe how the source catalogues 
in each band have been merged into a single catalogue, and 
how that catalogue has been used to obtain photometric red- 
shifts. These are then used in Section [4| to show which ob- 
jects are associated with the structure around the QSO and 
to the submm sources, and how we have obtained the dif- 
ferent physical parameters for each object. The evolutionary 
stage of our sources is discussed in Section [5] Finally, in Sec- 
tion [6| we summarize our results. We have assumed through- 
out this paper a Hubble constant Ho — 70 km/s/Mpc, and 
density parameters fini ~ 0.3 and SIa = 0.7. The spectral 
index a is defined as cc i^" . 



2 DATA 

The _R-band, K-h&nA, Spitzer and SCUBA data used here 



are the same ones discussed in Stevens et al. (2004 20101 
(henceforth SIO), but with a new reduction of the 7?-band 
data, all the astrometry tied to this band (see below) and 
fluxes for fainter Spitzer sources. The i, Z and J-bands and 
radio data are newly presented here. We show in Fig. [Tj the 
central area of the submm images along with contours and 
images of the radio data. The region to the NW around 
source 850_3 is shown in Fig. [5] 

We give in Table [l] the deboosted 850 ^m and 450 
fluxes from SIO, where we have symmetrized the error bars 
taking the worst of the upper and lower error bars. The 



source extraction method used in SIO (described in Scott 
et al. 2002 1 can miss real sources separated by less than a 



beam (14.2" for 850 ^m). This is probably the case for 450_3, 
which is blended together with 850_2 (see Fig. [l] and SIO): 
in 850 ^m these two sources are part of a ~N-S struture 
with 850_2/450_2 in the Southern tip and 450_3 close to the 
Northern tip. 

A summary of the origin of the observed optical-to- 
radio data used here is included in Table [TJ as well as some 
of the characteristics of the final images. The data come 
from several different telescopes and observatories. At all 
wavelengths shorter than radio, the data were taken as in- 
dividual images which were later combined using standard 
procedures (whose details may differ from data set to data 
set and are outlined below) . Since this involved combination 
of images, often with non-integer pixel shifts and individual 
re-scaling, we have used a simulation technique to take into 
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Figure 1. Finding cliarts of a 2.1' X 1.3' region around the submm structure associated with RX J0941 (N is up and E to the left). The 
region around 850_3 (to the NW of RX J0941) is shown sepparately in Fig. [s] Top left: 450 /.tm SNR image with 2 to 6(t contours (in 
steps of 0.5(7, dashed cyan) and 6cm 2 to 5ct contours (in steps of Ic, yellow), the numbers and crosses in cyan correspond to the 450 
sources in Table[l] Top right: 850 fim SNR image with 2 to 8a contours (in steps of O.So", dashed magenta) and 20cm 2 to 5a contours 
(in steps of la, yellow), the numbers and "x" in magenta correspond to the 850 fim sources in Table [T] (the cyan crosses to the 450 fim 
sources in that Table). Bottom left: 6cm image and contours (as above) with source numbers and crosses as above (note that sources 
1 and 2 are common to 450 ^m and 850 fiia). Bottom right: 20cm image and contours (as above) with source numbers and markers as 
above. In greyscale, magenta appears slightly darker than cyan and the steeper yellow contours appear white. 



Table 1. SCUBA source positions and fluxes (see Section [2]|. The first part of the table corresponds to 850 /.tm and the 
second to 450 fim. Our source number for the counterpart is given under # (see Table[2]and Section [4.2[ l, as well as the best 
fit photometric redshift (zp) and its la confidence interval interval. Star formation rates (SFR) and 8-1000/jm IR luminosities 
are from template SEDs (see Section |4.2[ l. Dust mass ranges and galaxy stellar masses are derived in Sections |4.4| and |4.5| 
respectively. Infrared luminosities, SFR, dust masses and galaxy stellar masses are obtained assuming z = 1.82 



Source 


R.A. 


Dec. 


Flux 


# 


Zp 


la 


-^^8-1000^111 


SFR 


Mdust 


logM.* 




(h m s) 


(°"') 


(nijy) 








(lOi^Lo) 


(Mo/yr) 


(lO^Afo) 


(logA/g) 


850_li 


09 41 44.63 


+38 54 38.88 


12.4+1.1 


39 






3.9 


4800" 


8-25 




850_22 


09 41 46.31 


+38 53 56.94 


5.7+2.1 


140 


1.85 


1.83-2.16 


1.3 


1700 


3-12 


11.6 


850-3 


09 41 42.62 


+38 55 36.81 


3.1+2.2 


501 


1.85 


1.83-2.20 


0.5 


650 


2-6 


11.5 


450_li 


09 41 44.57 


+38 54 39.46 


44.7+10.2 


39 






2.4 


2900 


5-22 




450-22 


09 41 46.31 


+38 53 56.22 


41.2+10.0 


140 


1.85 


1.83-2.16 


2.2 


2700 


4-23 


11.6 


450_3 


09 41 46.35 


+38 54 15.44 


33.5+8.8 


262 


2.78 


2.58-3.00 


1.8 


2200 


4-16= 


10.9^= 


450-4 


09 41 43.63 


+38 54 24.77 


28.8+7.8 


98 


1.84 


1.81-2.08 


1.4 


1700 


3-16 


11.5 



850-1 and 450_1 are the same physical source: the central Q SO. ^ 850-2 and 450-2 are the same physical source 
" Above maximum Ljr luminosity for |Chary fc Elbaz| l |200ll l 
^ Stellar masses are for the counterparts under column "#" 

" The best fit redshift is 2.78, very different from 1.82. For the higher redshift SFR=26OOM0/y, Lg-iooOMm = 2-1 X lOl^L©, 
log(Mt/MQ) = 11.3. Using the la uncertainty in the photometric redshift the uncertainty interval becomes 
Mdurt = (4-25) X IO^Mq 
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Table 2. Summary of the observations on the field around RX J0941 



Telescope 


InstmrQent 


uare 


Filter / channel 


T<-i-v-i-» 

lexp 


mo 


IcT detection level 


Seeing 












(Vega) 




(FWHM, ") 


WHT 


PFI 


28-May-2003 


R Harris 


2880 


32.35 


25.401 


0.993 


Gemini-North 


CMOS 


13-Feb-2005 


i 


3000 


34.34 


25.851 


0.743 


INT 


WFC 


03-Mar-2006 


Z 


72000 


29.49 


24.85^ 


1.275 


UKIRT 


UFTI 


ll-Feb-2004 


J 


7200 


30.18 


23.02^ 


0.741 


UKIRT 


UFTI 


24-May-2003 


K 


7700 


28.90 


21.87^ 


0.708 


Spitzer 


IRAC 


06-May-2005 


4.5 /im 


3000 




0.52 




Spitzer 


IRAC 


06-May-2005 


8 fim 


3000 




2.12 




Spitzer 


MIPS 


12-Apr-2005 


24 /xm 


336 




522 




VLA 


C band 


25-Mar-2008 


4860 MHz 


64000 




6.6^ 


_4 


GMRT 


L band 


21-Jan-2008 


1280 MHz 


25000 




29-^ 


_5 



1; in Vega magnitudes, from App. fusing 1.5 arcsec aperture and the "standard" aperture correction 
2: in ^Jy, see App . |A| 
^: in /iJy, see Sec. l2| 

beam size 3.3" X 2.2", position angle 57 deg 
^: beam size 3.6" X 3.5", position angle 73 deg 



account the effect of this procedure on the overall photomet- 



ric calibration (see Section 2.1 1 



The astrometry of the _R-band image was matched to 
APlvJ^ The astrometry of the iZJK and Spitzer images was 
then refered to the i?-band image. The astrometry of the 
SCUBA images in SIO was also tied to the i?-band image, 
allowing for a shift, asymmetric X and Y scales and rotation 
with respect to the "native" SCUBA astrometric system. 
The sources used as reference for these were 850_1, 850_2 and 
850_3 for the 850 /xm image and 450_1, 450_2 and 450_4 for 
the 450 /xm image. 450_3 was excluded because it is probably 
a blend of several sources (see above). 

The radio data have been reduced with AIPS follow- 
ing a standard procedure. For the VLA data we have used 
an intermediate weighting scheme between natural and uni- 
form (ROBUST=0 in AIPS). This gives a synthesised beam 
of 3.7" X 3.5" and a noise of 6.6 /ijy/beam. We have used 
the original astrometry for the radio images. The very good 
angular resolution and astrometry of the radio images allow 
relating better the submm detections to the corresponding 
optical counterparts. 

Finally, in Section |4.3| we have used the upper limit 
server FLIXp] (which in turn uses the procedure outlined in 
Carrera et al. 20071 to obtain la upper limits to the 0.5- 



2 keV (soft) and 2-12 keV (hard) X-ray countrates from 
the corresponding background maps on the positions of the 
SCUBA sources. These X-ray data are the same ones dis- 



cussed in Page et al. ( 2010 1, but the server uses the standard 



output from the XMM-Newton pipeline. 



2.1 Photometric calibration 



We took advantage of the existence of SDSS ( [Adelman 
McCarthy et al.||2007ll data in the RX J0941 field to tie 



1 'http: //www. ast . cam. ac.uk/~apmcat' 



http: //ledas-www. star . le . ac .uk/f lix/f lix.html 



[http:// 



www.sdss 



•°rg/f 



our optical photometry to the SDSS DR5 photometric cal- 
ibration. In summary, we obtained Vega magnitudes for 
the SDSS sources within the FOV of each of our images 
(using expressions from Jester et al. 2005 and from the 



SDSS and WFC survey web pages), using only point-like 
objects (type=6 in the SDSS source classification). We then 
matched the SDSS sources to the sources in our images, fit- 
ting an additive constant between the magbest SExtrac- 
TOR magnitudes of our sources (excluding saturated ones) 
and the Vega magnitudes, taking into account the statistical 
errors of both sets of magnitudes. These additive constants 
are the zero points given in Table [l] The error in the zero 
point is added in quadrature to the error in the magnitude 
of each source to get the final error in the magnitude of each 
source. 

Unfortunately, the FOV of the JK images was very 
small and there were no 2MASS sources within them to 
provide an independent check on the magnitudes. However, 
the sky conditions while our NIR images were taken were 
photometric and the templates fits (see Section [3| do not 
show any systematic trend for the NIR magnitudes. 

The Spitzer source lists in SIO had a high significance 
detection threshold, tailored to find counterparts to the 
submm sources. However, the deep RiJK images used here 
detected much fainter sources, which could be seen in the 
Spitzer images but were not present in those source lists. In 
order to have source lists which would include those fainter 
sources and to keep compatibility with the results of SIO we 
used S Extractor with lower significance to get aperture 
magnitudes (1.3, 1.7 and 3.2 arcsec radius for 4.5, 8 and 
24 ^m respectively, chosen to avoid source confusion) for 
the sources, and then obtained an empirical factor (taking 
into account the errors on both axis) which would match 
these SExTRACTOR fluxes with those of SIO for the com- 
mon sources. The SExtractor fluxes were then multiplied 
by these factors, taking also into account the uncertainties 
in the factors to calculate the flnal flux uncertainties. These 
flnal Spitzer fluxes are hence corrected to "inflnite" aper- 
ture. 
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Table 3. Positions and photometry (Vega magnitudes and fluxes) of tfie counterparts mentioned in the text 



M- 

it 


R A 


Dec 






7 

£j 


T 

o 












(h m s) 


(°"'') 












(MJy) 


(MJy) 




13 


09:41:45.84 


+38:53:53.54 


22.27+0.10 


21.42+0.07 


20.95+0.04 


19.87+0.18 


18.12+0.11 


37+2 


15+2 




16 


09:41:46.14 


+38:53:55.71 


23.36+0.17 


22.76+0.09 


22.82+0.17 


21.12+0.25 


19.07+0.13 








29 


09:41:46.17 


+38:54:17.99 


23.27+0.18 


22.84+0.10 


22.82+0.18 


21.50+0.33 


19.95+0.23 


13+1 






39 


09:41:44.64 


+38:54:39.57 


19.63+0.06 


18.73+0.07 


18.79+0.02 


17.70+0.17 


16.46+0.10 


715+36 


1661+87 


5443+360 


98 


09:41:43.56 


+38:54:24.33 


23.70+0.51 


23.00+0.18 


22.20+0.21 


20.59+0.30 


18.26+0.13 


58+3 


37+3 




135 


09:41:45.94 


+38:54:14.66 


24.56+0.52 


23.32+0.13 


22.84+0.17 


21.41+0.30 


19.51+0.16 


19+1 


13+2 




140 


09:41:46.31 


+38:53:56.55 


23.21+0.23 


22.70+0.12 


21.96+0.12 


20.07+0.20 


18.03+0.11 


85+4 


54+3 


699+69 


262 


09:41:46.35 


+38:54:13.94 


25.26+0.96 


25.88+1.09 




26.80+34.38 


21.13+0.54 


14+1 


14+2 




501 


09:41:42.65 


+38:55:36.37 


22.95+0.09 


22.86+0.08 


22.08+0.07 






69+4 


51+3 


345+57 


512 


09:41:42.53 


+38:55:36.85 


23.84+0.16 


22.67+0.08 


22.82+0.12 













Table 4. Positions and fluxes for the radio sources from VLA (4.86 GHz) and GMRT (1.28 GHz), as well 
as the spectral index a. The upper limits quoted are 3a-. 



Source 


R.A. 

(h m s) 


AR.A. 

(s) 


Dec. 

(°"') 


ADec. 

(") 


1.28 GHz flux 

(MJy) 


4.86 GHz flux 
(A'Jy) 


a 


850_1/450_1 

850_2/450_2 

850_3 

450_3 

450_4 


09 41 44.66 
09 41 46.32 
09 41 42.64 
09 41 46.39 
09 41 43.59 


0.01 
0.01 
0.03 
0.03 
0.01 


+38 54 39.9 
+38 53 56.7 
+38 55 36.8 
+38 54 13.8 
+38 54 24.3 


0.1 
0.1 
0.3 
0.3 

0.1 


650+29 
674+52 

<96 
<81 
189+29 


194+10 
313+11 
28.0+6.3 
28.5+6.1 
71.8+6 


-0.90+0.05 

'-'■'^'-0.06 

>-0.92 
>-0.78 

-0 y^+o-" 


Blob 


09 41 45.76 


0.03 


+38 54 07.8 


0.3 


800+100 


120+20 


_1 42+0-18 



3 PHOTOMETRIC REDSHIFTS 

Once -we have the final images in each band, we have run 
S Extractor on each one of them with fairly low signif- 
icance requirements (3 or more pixels 3 or more a above 
the background) to obtain comprehensive source lists in 
each band. These source lists have been restricted to the 
FOV of the JK images (roughly 1.5 arcmin squared), to en- 
sure the maximum wavelength coverage for the photometric 
redshifts. If a source was detected in either of the Spitzer- 
IRAC bands, we have included the corresponding measure- 
ments, otherwise, we have just used the 5 optical-NIR bands 
RiZJK to assign the photometric redshifts. The individual 
RiZJK and IRAC source lists have been merged, consider- 
ing two sources to be the same if they are closer than 1 arcsec 
to each other. The positions, magnitudes and fluxes of the 
counterparts mentioned below are given in Table [2] 

We have a total of 239 unique sources detected in at 
least one of our 8 bands. Photometric redshift fits could not 
be performed on 14 of them, mostly because at least four 
of the five bands were outside the FOV and, in two cases, 
because the two measurements which were within the FOV 
were upper limits. 

To obtain photometric redshifts we have used hyperz 



(version 11, Bolzonella, Miralles, & Pello 2000 1. We have not 
allowed for intrinsic reddening but we have included Galac- 
tic de-reddening {E{B — V) = 0.015). We have allowed a 
redshift interval between and 3. Many current photomet- 
ric redshift codes (e.g. |Bemtez||2000[ ) take into account the 
expected apparent magnitude distribution of the objects (us- 
ing galaxy luminosity functions), reducing the likelihood of 
redshifts that would result in absolute magnitudes too bright 
or too dim for the observed magnitudes, hyperz does not 



provide this facility, so we have instead checked "a poste- 
riori" the absolute B-band magnitude of the objects. Fol- 
lowing [Rwan^Robinso^^eFal^ henceforth MRR08, 
we have adopted redshift-dependent lower and upper lim- 
its to the absolute blue magnitude of max(-19.5,-17-z) and 
max(-25,-22.5-0), respectively. 



We have tried several sets of templates (Bruzual & 
Chariot |[2003l |Rowan-Robinson at ar] |2008) and settled for 



the "new" galaxy templates described in MRR08. There are 
7 different templates (starburst SB, youngE, E, Sab, Sbc, 
Scd, Sdm). The last 5 templates define a sequence of in- 
creasingly later galaxies, and template youngE is a "young" 
Elliptical with an age of IGyr, to allow for the limited time 
for evolution available at higher redshifts. These templates 
do not include dust emission, so they are appropriate until 
about 3 fim rest-frame, including our observed RiZJK and 
IRAC bands at the relevant redshifts. When comparing the 
full SED we will complement those with other SEDs with a 



wider coverage (see Section 4.2 1 



For a good match of the spectral shape of the templates 
to the observed photometric points it is essential that pho- 
tometry of each source corresponds to the same aperture. 
To minimise the mutual photometric contamination by close 
sources in the optical-NIR bands, we have used an adaptive 
aperture size with a maximum radius of 1.5" (see App. [a|. 

Out of the 225 sources in which photometric redshift fits 
could be performed, if we select those with a 99% redshift 
interval that includes z = 1.82 and that is narrower than 
0.6 (full span from lower to upper limits, i.e., relatively well 
constrained) , we get 4 sources including all the counterparts 
to the SCUBA sources (except for 450_3). 



Using the Bruzual & Chariot ( 2003 1 templates instead 
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(which also include blue templates with recent bursts of star 
formation) the best fit photometric redshifts for the coun- 
terparts to 850-2, 850-3 and 450-4 are the same as with the 
MRR08 templates, and these 3 sources are also among the 
8 with "narrow" 99% confidence interval including z — 1.82 



with the Bruzual & Chariot ( 2003 1 templates 



We conclude that our photometric redshifts are robust 
against using different sets of templates. In particular, the 
counterparts to the SCUBA sources (see below) always ap- 
pear amongst the sources with higher quality fits compatible 
with the redshift of the central QSO. 



4 RESULTS 
4.1 Radio data 

The list of the radio sources detected can be found in Ta- 
ble [3] The very good alignment of the radio sources (with 
their original astrometry) to the submm sources (aligned in- 
dependently to the J?-band image) can be seen in Figs.[T]and 
[5] All submm sources are detected with high SNR> 4 in the 
radio images (except for 450-3, which also has the optical- 
NIR counterpart less likely to be associated with the QSO, 
see below). This increases our confidence both in the exis- 
tence of the submm sources and on their association with 
particular optical-NIR sources, since the angular resolution 
of the radio images allows very accurate positioning of the 
long wavelength sources. 

It is worth noting the existence of a "blob" of diffuse ra- 
dio emission to the SW of 450-3 (09:41:45.694 -^38:54:07.57), 
clearly visible in both the VLA and GMRT data with a di- 
ameter of about 3 arcsec. It appears to be connected to 450_3 
by a bridge of ~ 2a emission in the former image. Its spec- 
tral index is a = —1.65 ± 0.20 which would advocate for it 
to be a lobe of radio emission. This area of diffuse emission 
is also clearly devoid of optical and Spitzer counterparts, 
which would be consistent with the lobe hypothesis. There 
is no evidence of a counter-lobe to the NE of 450_3. If it is a 
lobe, it could be associated instead with source 850-2/450-2, 
which is much brighter in radio, but no bridge of emission 
is seen between the latter source and the "blob" . 



4.2 Association to submm sources, SEDs and SFR 

We have plotted finding charts in i, K, Spitzer, SCUBA, 
VLA and GMRT for the 5 distinct sources from SIO (see 
Figures |3] to [7|. The best fit (to RiZJK, 4.5 and 8 fim) 
photometric redshift templates are shown as red lines in Fig. 
[2] as well as the observed wide-band SEDs with a number of 
templates. All the redshift confidence intervals given in the 
rest of this Section are la. 

The templates from MRR08 do not include dust emis- 
sion, so in order to model adequately the longer wavelengths 
we have used the starburst models of |Chary fc Elba^ ( |2001[ ) 
(henceforth CEOl), which have been obtained from observa- 
tions at several NIR-FIR wavelengths requiring the relative 
fiuxes in different bands to be consistent with the observed 
IR luminosity functions. One feature of these models which 
is worth noting is that their shape is luminosity-dependent 
and hence, given a redshift and the flux in one band, the SED 
is fully defined. The authors have provided IDL routines to 



compute their templates at a given redshift and through 
a number of instrumental bandpasses. We have added the 
SCUBA 450 and 850 /im sensitivity curves to the set of avail- 
able bandpasses. The absence of a good correlation between 
the B-band and IR luminosities implies that the optical/NIR 
part of their SEDs is uncertain. This is not a major problem 
here since we are only using the CEOl SEDs for the longer 
wavelengths. Instead of performing a formal fit to the 
data, we have computed the CEOl templates corresponding 
to each of the observed Spitzer and SCUBA fluxes at the 
redshift of each source. 

In the case of 850_1/450_1 we expect to have a strong 
AGN component, specially in the rest-frame optical/NIR 
part of the spectrum. We have modelled it using templates 
from MRR08. We have approximated their RRlv2 AGN 
template with a broken power-law (which is a good approx- 
imation to the continuum shape), redshifting it to 2: = 1.82 
(the r edshift of RX J0 941). Within the Unifled Model for 
AGN ( Antonucci|1993 1 part of this direct emission is inter- 
cepted by a dusty torus, which "reprocess" it and emits it 
at mid and far-infrared wavelengths. To model this emission 
we have used the dust torus template from MRR08. 

The SFR from the CEOl templates normalized to 
450/850 fim are given in Table [1] The SF R have been calcu- 
lated using the following expression from Kennicutt ( 1998 1 



SFR{Mq/j) = 1.7217 X 10-'°Lfir(Lq), 



(1) 



where L(8 - 1000/im) (Lir) and L(40 - 500^im) (Lfir) 
have been estimated from the 12, 24, 60 and 100 /xm 
fluxes (obtained from the CEIO templates normalized to 
the SCUBA fiuxes) using the formula described in Sanders 
& Mirabell (|1996|| with C = 1.4. Most of our submm 



sources would be classified as Hyper Luminous IR Galax- 
ies using the criteria in Rowan- Robinson ( 2000 1 . Since 



the CEOl templates are simply rescaled to the observed 
SCUBA fiuxes (F), the uncertainties can be estimated via 
ASFR/SFR=AF/F. 

We now discuss the counterparts to the submm sources 
and the SEDs individually: 

• 850-1/450-1: The sharp radio contours prove that this 
source is indeed the central QSO (RX J0941). It is asso- 
ciated with source 39 (242/180 in SIO). The fuU SED for 
this object is shown in Fig. [2] (top left panel). The dashed 
red line represents the CEOl template normalized to the ob- 
served 450 /im fiux at 2 = 1.82. It clearly reproduces well the 
observed emision at 450 /im and red wards, but it falls well 
short of the observed values bluewards of that wavelength. 
We have fitted the RiZJK fiuxes with our broken power-law 
model for the "direct" AGN emission from MRR08 (green 
dashed line), which more or less follow the blue spectrum in 
the observed optical/NIR part of the SED, but again fails to 
account for the observed MIR emission. We have added the 
starburst emission to the direct AGN emission, subtracted 
it from the 8 /im point and rescaled the dust torus emission 
from MRR08 to this value (blue dashed line). The sum of 
these three components is shown as the black solid line in 
the Figure. It clearly follows the overall shape of the SED, 
which cannot be reproduced by any two of the models alone. 
The submm emission is clearly dominated by the starburst, 
with the total AGN contribution being at least two orders of 



magnitude smaller, as originally stated in Page et al. (2001 1 
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Figure 2. Spectral Energy Distribution in F^, in the observer's frame for 850_l/450-l (RX J0941, top left), 850_2/450_2 (top right), 
850_3 (middle left), 450_4 (middle right), and 450_3 (with templates at the best photometric redshift -bottom left- and fixed at z = 1.82 
-bottom right-). Observed points (RiZJK, 4.5, 8, 24, 450 and 850 fj-ra, and 6 and 20 cm) are shown as black dots and upper limits as 
down-pointing arrows. The meaning of each line is summarized in each panel and explained fully in Section |4.2| 
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Figure 3. Finding charts for 850_1/450_1. From top-left to bottom-right: i (the size of the red box is 10" X 10"), K, 4.5 fim (dashed 
cyan circle is the positions of the 4.5 lira sources with radius equal to the flux extraction radius), 8 /im (dashed cyan circles mark 8 ^m 
sources and extraction radii), 24 fim (dashed cyan circles mark 24 /im sources and extraction radii), 450 /im (with cyan contours starting 
at 2(T and increasing in steps of 0.5), 850 fim (with magenta contours, starting at 2cr and increasing in steps of 0.5), 4.86 GHz (with 
yellow contours starting at 2a and in intervals of 1), 1.29 GHz (with yellow contours, same as previous). All finding charts have the 
same FOV and the solid circles indicate the optical/NIR source position and extraction radius (green corresponds to 1.5", red to smaller 
radii), only the sources mentioned in the text have been labelled. The labels have been omitted in the SCUBA and radio images for 
clarity, but we have labelled and marked instead with cyan "+" our positions for the 450 ^m and with magenta "x" the positions of the 
850 /im sources. 



There is very little room for an additional radio contribution 
from the AGN in this radio-quiet QSO. 

• 850_2/450_2: The radio data show a strong pointlike 
source, coincident with our source 140 (191/155 in SIO). Its 
photometric redshift is Zp = 1.85 with a redshift interval 
spanning 1.83-2.16 with spectral type Sab. Our source 16 is 
about 2.2 arcsec away from 140, and could be contributing 
to its 24 /im flux. Its photometric redshift is 3.0 (2.9-3.0, 
Scd). Source 13 is brighter and coincides with elongations in 
the contours at 450 /im and 850 /im. Its photometric redshift 
is Zp ~ 1.19 ± 0.03 (Scd), and hence is probably unrelated 
to the central QSO. In Fig. [5] (top right panel) the dashed 
and dotted lines show the CEOl templates normalized to 450 
and 850 /im (respectively). These lines straddle very closely 
all bluer and redder wavelengths. 

• 850_3: The NE-SW-elongated 2 and 2.5a 850 /im con- 



tours are aligned with the VLA 2a contour. However, there 
is no sign of it in the GMRT image. This peak would favour 
our source 501 as the prefered counterpart (zp = 1.85, z in 
1.83-2.20, Sab, 339/241 in SIO), ahhough this optical/NIR 
neighborhood is rather crowded: sources 502, 511 and 512 
are also close. The Spitzer emission appears pointlike favour- 
ing either 501 or 512. Unfortunately, the absence of J and K 
data for this area does not allow to put strong constraints 
on the photometric redshift range for 512 {zp ~ 0.7, SB). 
In Fig. [2] (middle left panel) the CEOl template normalized 
to the 850 /im flux (blue dotted line) reproduces reason- 
ably well all spectral points redwards of 8 fim, but falls well 
short at bluer wavelengths, although the spectral shape is 
similar. Normalizing the CEOl template to the 4.5 fim flux 
is a bit too high at 24 fim and 850 /im, and too high at radio 
wavelengths. 
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Figure 4. Finding charts for 850_2/450_2. Same layout as in Fig.js] 



• 450_3: The 450 and 850 fim contours are complex (this 
source might also be detected in 850 /im, see Section [2]| The 
> 2(7 VLA contour includes the 450 /xm peak of 450_3, co- 
inciding with our source 262 (206/160 in SIO), and hence is 
our prefered counterpart. Object 262 comes originally from 
a 4.5 detection and only has upper limits in iZJ, hence 
its redshift is poorly constrained, seeming to prefer Zp ~ 2.8 
(youngE, z > 2.6). Its absolute B magnitude would be too 
faint for both z = 1.82 and z = 2.8, suggesting an even 
higher redshift. This source seems to be connected by a 2a 
bridge of radio emission to the radio "blob" discussed in 
Section [4.1| Nearby sources with some evidence for Spitzer 
emission are 135 (205/159 in SIO, Zp = 1.85,1.8-2.2, Scd) and 
29 (to the NW and along a N-S elongation in the 850 /im 
contours, Zp — 2.7, 2.6-2.9, SB). Object 135 coincides with 
the main E-W asymmetry of the 450 fim emission and could 
be related to the central QSO, although the absence of a ra- 
dio detection does not support this idea. Since the redshift of 
262 is poorly constrained we show the CEOl SED normalized 
to the IRAC 4.5 fim and SCUBA points both at z = 2.78 
(bottom left) and at the redshift of the QSO z = 1.82 (bot- 
tom right). The CEOl template normalized to the 450 /im 
point is well above the radio and the 4.5 and 8 /um points in 



both cases, while when it is normalized to the 4.5 fim point 
it matches better the IRAC and the radio points. This may 
mean that there might be some contribution to the SCUBA 
flux from our source 135, in agreement with the elongated 
contours at 450 fim. The 24 ^m source position is interme- 
diate between 262, 135 and 29, so the flux in that range also 
probably includes contribution from them, explaining why 
it lies above the CEOl template at z = 2.78. 

• 450_4: Coincides with a well defined pointlike IRAC, 
VLA and GMRT source, very close to our source 98 (the 
unlabeled source to the W of the centre of the finding chart 
for 450_4 in SIO). We have found Zp = 1.85 (1.8-2.1, Sbc). Its 
SED is plotted in Fig. [2] (middle right panel), showing that 
the CEOl models at that redshift neatly follow the shape 
of the Spitzer and SCUBA points (only the one normalized 
to the SCUBA point is shown, the models normalized to 
the Spitzer points are indistinguishable in this scale). The 
only clearly discrepant range is radio, which falls a factor of 
about 3 below the CEOl template. 

In summary, the radio positions allow assigning un- 
ambiguous counterparts to all of the submm sources. 
Four out of five submm sources (850_1/450_1 -RX J0941-, 
850_2/450_2, 450_4 and 850_3) have counterparts with red- 




shifts compatible with 1.82 within ~ la. The three latter 
ones have very similar spectral shapes, compatible with an 
existing relatively old stellar population with ongoing star 
formation. The last submm source (450_3) appears to be a 
radio emitting object at higher redshift, although there is 
a nearby counterparts within the "right" redshift interval 
and very similar properties to the other counterparts to the 
submm sources (135). From the number counts in Coppin 



|et al.| p006| ) we would expect to have about one submm 
source unrelated to the QSO in our field, so it is not sur- 
prising that 450-3 turns out to be a "background" object, 
considering also its radio faintness. 



4.3 Upper limits on AGN luminosities and 
estimating black hole masses 

850-1/450-1 is the central QSO and hence it obviously har- 
bours a powerful AGN. In this Section we will concentrate 
on the rest of our submm sources. The "blob" of radio emis- 
sion close to 450-3 might be interpreted as evidence for the 
presence of an AGN in this object. This possibility cannot 
be rejected using our current data due to its quality at the 



bluer wavelengths, where the AGN would be expected to 
dominate. 

The optical-to-radio SEDs of the other three counter- 
parts to the SCUBA sources have very limited room left for 
an additional AGN emission component. We will derive be- 
low limits to the luminosity of any putative AGN and the 
mass of the associated SMBH, under the assumption that 
the latter is growing through accretion and radiation is emit- 
ted in the process. Our galaxies could harbour SMBH more 
massive than deduced below if they were dormant. 

We have taken advantage of our XMM-Newton obser- 



vation of this field (Page et al. 20101 to set limits on the 
intrinsic luminosity of any such AGN (see Sectio n [2|. Of 
course, the central QSO is actually detected (see [Page et| 
al.|2010| . Three other SCUBA sources have upper limits on 



the pn camera (850-3 is in a gap between two chips in that 
camera), with similar values of <0. 00034 (<0. 000455) cts/s 
in the 0.5-2 keV -soft- (2-10 keV -hard-) band (using a sim- 
ple power-law with photon index F = 1.9 and local Galac- 
tic absorption, these values correspond to fluxes of about 
< 6x lO^^'' and < 4x 10"^^ erg cm s , respectively). We 
have taken these countrates and estimated what would be 
the intrinsic luminosity of a 2 = 1.82 AGN to produce those 
countrates if viewed through different amounts of intrinsic 
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absorption, parameterized as the Hydrogen column density 
Nh- We have used Nh=0, and log(iVH/cm"^) = 21.5-25.5 
in steps of 1 dex. We adopt a simply parametrised AGN 



spectral model following Gilli, Comastri, & Hasinger (20071, 
i.e. an exponentially cut-off absorbed power law of photon 
index F — 1.9, a small scattered component, a reflection 
component, and an Fe K emission line. We include a correc- 
tion for Compton scattering of the continuum at the highest 
column densities. We have used the program xspec and an 
on-axis response matrix and ancillary file for pn to convert 
from countrates in those bands to luminosities in the rest- 
frame 2-10 keV band ix,2-iokcV- The results are shown in 
Fig-i 

Supernovae and binary stars associated with star for- 
mation produce X-rays, which would further reduce the al- 
lowance for any AGN emission. [Ranalli et al.| ( |2003[ ) found 
a tight correlation between X-ray luminosity and FIR lumi- 
nosity in local starburst galaxies which, using Eq. 1 above 
( |Kennicutt||1998[ ), can be expressed as SFR(MQ/y) = 2 x 
10~*"Lx,2-iokev(erg/s) (their equation 15). For the typical 
values of our galaxies of SFR~ 2OOOM0 /yr this corresponds 
to I/x,2-iokoV ~ 10*^ erg/s, comparable with the above up- 
per limits for the soft band and moderate absorption. 

In addition, we can set very conservative upper limits to 



the presence of an obscured AGN from the observed 24 fim 
flux. We have used as fiducial values those of sources 850_2 
and 850_3 (~700 and ~ 350/iJy, respectively) , because 450_4 
is not detected in that band and 450_3 is probably a back- 
ground source and its emission in that band is probably a 
blend of several sources. Assuming that the emission from 
dust heated by the AGN at 24 /xm is isotropic and that the 
full observed fluxes at that wavelength have this origin, we 
can use an AGN SED (the [Richards et"aL][2006| RQ QSO 
template redshifted to z = 1.82) to find out what hard X- 
ray fluxes those 24 nm fluxes would correspond to. The cor- 
responding values are log(Lx,2-iokev/ergs~^) = 44.3 and 
44.0, respectively. Since the starburst SED template leaves 
very limited room for an additional AGN contribution (for 
instance, for 850_2 the difference between the 24 ^m flux 
of CEOl template normalized to 450 /im and the actual ob- 
served flux at that wavelength is about 200 ^Jy, a factor of 
about two below our lowest fiducial value), this is an abso- 
lute ceiling to the maximum luminosity of a putative AGN 
in the innermost regions of our submm sources. 

We conclude that AGN with intrinsic 2-10 keV lu- 
minosity ~ 10*^ erg/s could be present in our SCUBA 
sources, if they are behind modest amounts of absorp- 
tion (log(7VH/cm~^) ^ 22 — 23), while objects with al- 
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most QSO luminosities (~ 10*'' erg/s) would have to be 
hiding behind much more substantial amounts of mate- 
rial (log(iVH/cm~^) > 24). More luminous AGN are not 
compatible with our current data. An X-ray luminosity of 
~ 10** erg/s is compatible with the observed hard X-ray 
luminosities of SMGs at similar X-ray fluxes ( [Alexander et| 
aLl|2005| 



We now translate these limits on X-ray luminosity into 
a limit on the putative SMBH mass. The first step is to esti- 
mate the total bolometric luminosity of the AGN, which we 
can write as Lagn = Kix,2-iokoV5 where k is the bolometric 
correction. The maximum theoretical luminosity emitted by 
material falling into a black hole of mass M, is called Ed- 
dington luminosity LEdd ( |Eddington|[T9T3l [Rees|p84| and 
can be expressed as 



log 



Lr. 



4.53 + log 



M. 



(2) 



Of course, a source can be emitting below this limit. The 
Eddington ratio rj is defined as = Lagn jl^Edd, and hence 



log 



M. 



loe 



X,2-10kcV 



- 4.53 4- log(K/r7). 



(3) 



Replacing the upper limit found above (Lx < 10** erg/s= 
10'° *^ Lq) we finally get 



log( ^ ) <5.89 + log(K/»7). 



(4) 



Alexander et al. ( 2008 1 study the black hole masses. 



Eddington ratios and X-ray luminosities of obscured SMGs 
at z ~ 2, finding that for objects with Lx ~ 10** erg/s, 
log(M./M0) ~ 7.8 with t] ~ 0.2. They also find similar 
values of r\ for nearby obscured ULIRGs. Finally, the gas- 
rich environment in which the SMBH are likely to reside 
and the large fraction of SMGs that host AGN activity (see 
Alexander et al.]|2005| [Borys et al.||2005l but see also |Laird| 
et al.|20To I suggest that the accretion is likely to be reason- 
ably efficient, so we adopt a fiducial value of = 0.2. Using 
the luminosity-dependent relation between X-ray luminos- 



ity and bolometric luminosity of Marconi et al. ( 2004 1, gives 
K = 33, very similar to the value of k = 35 which would 
be ded uced from the QSO SEP of [Elvis et al.| |T994l), as 



used by Alexander et al. ( 2008 1. Vasudevan & Fabian ( 2007 \ 



find K = 15 — 25 for AGN emitting at low Eddington ratio 
(77 < 0.1) and «: = 40 — 70 for higher Eddington ratios. We 
will assume k = 35, getting 
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Figure 8. Upper limits to the 2-10 kcV intrinsic luminosity of 
an AGN from the upper limits to the soft band (black solid line) 
and hard band (red dashed line) count rates from XMM-Newton 
observations versus the intrinsic column density. The value for 
Nh = has been plotted at \og{NH/cm~^) = 20. The blue 
dotted lines mark the limits from the observed 24 /im fluxes for 
sources 850_2 and 850_3, see text for an explanation. 



log 



M. 



< 8.1 + log — 



(5) 



Our logarithmic upper limit on the SMBH masses can be as 
high as 8.4 {-q — 0.2, k = 70) or as low as 7.4 (rj — — 
35) for efficient accretion. For inefficient accretion (assuming 
K = 15) we get 8.0 for r] = 0.1 and 9.0 for r/ = 0.01. 



4.4 Dust masses 

We have estimated dust masses from the observed SCUBA 



fluxes using Eq. 5 in Martmez-Sansigre et al. ( 2009 1 



(1 + z)k(i/(1 + z))B,{T, v(\ + z)) ' 



where 



k{i'(1 + z)) = kq X 



iy(l + z) 



(6) 



(7) 



and Bv{T,v) is the usual blackbody emissivity (e.g. Eq. 3 
in Martmez-Sansigre et al.|2 009). 

For Fv we have used the monochromatic fluxes at 
A=450 and 850 ^m from Table [l] In order to estimate a 
range of dust masses for each object, we have used the eight 
possible comb inations of (kq = 0.04 m'^/kg, Aq =1.2 mm) 



(Beelen et al., 2006) and (kq = 2.64 mVkg, Ao =125 ^.m) 



I Dunne et al.j2003[ ) with (r=47K, /3 = 1.6) (representative 
of unobscured quasars at high redshift [Beelen et al.||2006[) 
and (r=35K, /3=1.5) (representative SMGs; see Kovacs et 
al. 2006) and with the minimum and the maximum of the 
Icr uncertainty interval in their photometric redshifts, taking 
the minimum and the maximum values of those combina- 
tions (which always corresponded to the minimum redshift, 
higher kq, higher T and maximum redshift, higher kq, lower 
T, respectively). We have made two exceptions: 850_1/450_1 



(the QSO) for which we have fixed z = 1.82, and 450_3 for 
which we give both the intervals fixing z = 1.82 (main body 
of the table) and allowing the full photometric redshift un- 
certainty interval (footnote). 

An additional source of uncertainty comes from the er- 
rors on the submm fiuxes, which range from about 10% to 
about 40% (with one case reaching 70%). These fractional 
errors are generally smaller than the uncertainty ranges de- 
rived above, and therefore we take the latter to represent a 
good estimate of the order of magnitude of the dust masses 
present in our sources. 

For the sources detected at both submm wavelengths, 
we find that the dust mass estimates derived from 450 and 
850 /im are compatible within the uncertainties. All sources 
have dust masses in the range 2 x 10* — 3 x lO" solar masses. 

We can also estimate the gas mass present in the central 
regions of those galaxies from the dust masses, assuming a 



gas-to-dust ratio of 54, deduced by [Kovacs et al. ( 2006 1 for 
z = 1 - 3 SMG (with an uncertainty of about 20%). This 



value is similar to the one obtained by Seaquist et al. ( 2004 \ 



for the central regions of nearby SCUBA galaxies. For our 
typical dust mass of ~ 10^ Mq this corresponds to a gas 
mass of ~ 5 X 10^°Mq. 



4.5 Galaxy stellar masses 

We have estimated stellar masses from the absolute ii'-band 
magnitudes Mk using the expression 



log 



M, 



Mk - 3.28 

^2^5 



log(3.2). 



(8) 



which we have taken from Borys et al. ( 2005 I assuming a 



light-to-mass ratio Lk/M = 3.2 x Lk,q/Mq, which is in- 
termediate between the average values from their fits to a 
burst and a continuous star-forming model. 

We have obtained Mk at the best photometric redshift 
integrating the best fit template using hyperz. These esti- 
mates should be reasonably accurate for the sources with 
IRAC detections up to z ~ 2.6, since the best fit template 
would then be "anchored" in a photometric point on a neigh- 
bouring rest frame wavelength. 

No estimate of the uncertainty in the normalization of 
the best fit template is returned by hyperz. Our first ap- 
proximation to the uncertainty in our derived stellar mass 
comes from the uncertainty in the redshift: 



Alog(M.) 



AMr 
2.5 



2Alog(di(z)), 



(9) 



since Mk = — 5 log(dL (2)) (distance modulus). We have 
estimated the uncertainty in the luminosity distance as half 
the difference between the luminosity distances at the lower 
and upper limits of the la redshift interval. 

Since essentially the derived stellar mass is proportional 
to the flux F at the observed wavelength corresponding to 
the rest-frame K-hand (between 4.5 and 8 /im at z = 1.82), 
a second contribution to its uncertainty would be: 



Alog(M.)' = Alogi^ : 



1 AF 



(10) 



In 10 F 

This corresponds to about 0.03, 0.04 and 0.07 for 450_4, 
850_2 and 850_3, for total errors of 0.09, 0.10 and 0.13, re- 
spectively. 
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A further source of uncertainty is the assumed light- 
to-mass ratio. MRR08 use a different approach to estimate 
the stellar masses from their best fit templates 3.6 /im lu- 
minosities, with a light-to-mass ratio that depends on the 
galaxy class and on the cosmic time. We have also esti- 
mated the stellar masses of our 3 counterparts following this 
recipe. We obtain L/M = 0.05 - 0.06 and log {M,/Mq) = 
11.85,12.06,12 for 450_4, 850_2 and 850_3, respectively. 
These stellar masses are between 2 and 3 times larger than 
those in Table [T] the discrepancies being larger than the 
total error budget calculated above. 

In what follows, we will use for our sources a fiducial 
value log(M*/M0) = 11.5 ± 0.2 which takes into account 
the different stellar mass values for the different sources and 
the redshift and photometric uncertainties discussed above. 
Using the MRR08 recipe would produce values higher by 
0.3-0.5 dex. 

These stellar masses are similar to those obtained 
for other samples of SMGs: e.g., Borys et al. (20051 
find log(M,/M0) = 11.4 ± 0.4 in the CDF-N/GOODS- 



N region, and Michalowski, Hjorth & Watson ( 2010 1 find 
log(M,/M0) = 11-12 with a median of 11.7, despite us- 
ing very different SEDs and methods for deriving the stel- 



lar mass (in the second case). As pointed out by Borys et 
al. (20051 those stellar masses are about 10 times larger 



than those of typical UV-selected star-forming galaxies at 



similar redshifts (Shapley et al. 20051. They are compat- 
ible with Schechter stellar mass for the GOODS-MUSIC 
galaxies at z ~ 1.8 (log(M,/M0) 



11.3 ± 0.15 Fontana 



et al. 20061. They also similar to the Schechter mass of 



3.6-4.5 ^m-selected galaxies in the HDF-N, CDF-S and 
the Lockman Hole in the 1.6 < 2 5C 2.0 redshift interval 
(log(M,/M0) = 11.40 ± 0.18 [Perez-Gonzalez et al.|[2008|. 



Their < 2 < 0.2 galaxies have slightly smaller masses 
(log(M*/M0) = 11.16 ± 0.25), although compatible within 
la. 

Although the above masses are stellar masses for the 
whole galaxy, they are probably dominated by the spheroid 
(as discussed by Borys et al.||2005 Alexander et al.||2008 l, 
since high-angular resolution H-band imaging of SMGs sug- 
gests that the stellar structure of most SMGs is best de- 
scribed by a spheroid/elliptical galaxy light distribution 
( Swinbank et al.|2010[ ). 



5 EVOLUTIONARY STAGE OF THE 
COUNTERPARTS 

In this section we use our estimates of the SFR and the stel- 
lar and dust mass for the counterparts to the submm sources 
to investigate their current state, their plausible future de- 
velopment and their relation to any putative AGN. 

An immediate conclusion is that our submm galaxies 
are very mature; they have stellar masses typical of MIR 
and submm-selected galaxies (~ 3 x lO^^Af©), and their gas 
masses (~ 5 X lO'^Afe) are around an order of magnitude 
smaller than our estimates of the stellar mass, indicating 
that >90% of the maximum possible stellar mass is already 
in place. 



Interestingly, Onodera et al. (20101 find a galaxy at a 



very similar redshift (2 = 1.82) with a very similar stel- 
lar mass, which has properties fully consistent with those 



expected for passively evolving progenitors of present day 
giant ellipticals. When the vigorous star formation taking 
place in our SMGs is over, they might well look like that 
object. 

A common indicator of the strength of the star forma- 
tion in galaxies is the Specific Star Formation Rate (SSFR) 
defined as the ratio of the SFR to the stellar mass. It 
is generally found that high mass galaxies tend to have 
lower SSFR at 2 < 2 (see e.g. Feulner et al. 2005 Perez- 



Gonzalez et al.|2008 Ilbert et al.|2010 1 , which is interpreted 



as a sign that star formation does not significantly change 
their stellar mass over this redshift interval. Our SMGs 
have SSFR~ 3 - lOGyr'V These values are about an or- 
der of magnitude higher than those of MIR-selected galax- 
ies with similar redshifts (1.6 < 2 ^ 2) and stellar masses 
(log(MH,/M0) > 11) in [Perez-GonzAlez et al] (|2008[). Our 



SMGs are amongst the most actively star forming objects 
of their time. 



Ilbert et al. (20101 find that, for a fixed range in SSFR, 



the fraction of high mass objects among their starforming 
galaxies drops with redshift between 2 = 2 and 2 = 1, 
while that of lower mass objects is approximately constant 
down to 2 ~ 1. This implies that the low mass star- forming 
galaxies are able to maintain a high SSFR, while the mas- 
sive galaxies evolve rapidly into systems with a lower SSFR 
(since the mass cannot go down, this means necessarily lower 
SFR). This is the expected evolution of our SMGs as de- 
duced above from an independent argument, and follows a 



clear "downsizing" pattern (e.g. [Gowie et al. 
Gonzalez et~al]|2008[ ) 



1996 Perez- 



Assuming that our SMGs host a SMBH which is grow- 
ing and emiting radiation, in Section 4.3 (equation [5| 
we deduced a very conservative upper limit to the mass 
of the SMBH in AGN in the centre of our galaxies of 
\og{M,/MQ) < 8.1. Hence, their black hole mass-to-galaxy 
mass ratio is log(M,/Mt) < —3.4. Applying the rela- 
tion between those quantities found by [Marconi fc Huiit] 



( 2003[ ) for local galaxies we get log(M./M.) = -2.6 (for 
log(M«/M0) = 11.5), so at face value our sources seem to 
have black hole masses at least about a factor of 6 below lo- 
cal galaxies of the same mass. Other studies have also found 
lower SMBH-to-galaxy mass values in 2 ~ 2 SMG (e .g. [Bo- 
rys et al.|[2005l [Alexander et al.[[2008[ [Coppin et aL][2009[ ), 



which are commonly interpreted as the black hole growth 
lagging behind the galaxy growth, since galaxies appear to 
be essentially fully grown while AGN are "under size". 

There are considerable uncertainties in such estimates 



(as discussed at some length by e.g. Borys et al. 2005 



Alexander et al. 20081. We include 0.2 dex of uncertainty 



in the stellar masses of our SMGs and the upper limit to 
the SMBH mass ranges between about 7.4 and 8.4 (in log- 
arithmic solar masses) for different Eddington ratios and 
bolometric corrections (see Sections [4.3[ and [4.5[ ). Taking 
the extreme values for both masses, we get upper limits to 
log(M,/M,) between -2.9 and -4.3, while the value for lo- 
cal galaxies is log(M./Mt)iocai = —2.6 ± 0.3 (using the rms 
intrinsic dispersion quoted in [Marconi fc Hunt[|2003l using 
instead [Haring fc Rix|2004| we would get -2.5 ± 0.3). 

Hence the upper limit to the ratio we derived above 
between the local and our SMBH-to-galaxy mass is 
log(M./M.)iocai/log(M./A/,) > 6, with an attached "un- 
certainty interval" ~ 1 — 100. The AGN in the centres 
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of our galaxies would have to be heavily obscured and 
just at the limit of detection for that ratio to be close 
to unity. The X-ray obscured AGN in SMGs studied by 



of radio emission. The photometric redshifts appear to be 



Alexander et al. ( 2005) are indeed heavily obscured 
have log(A'^ff /cm~ ) > 23) but their average luminosity is 
ix,2-iokeV ~ 5 X 10*^ erg/s, leading to an additional factor 
of 2 on that ratio. Finally, another factor of 2 or 3 would be 
necessary if we used the MRR08 method for estimating the 
stellar masses. 

In summary, our SMGs appear to be mature 
log(M*/M0) = 11.5 ± 0.2 and with limited scope of sig- 
nificant further increase in mass (remaining gas mass ~ 
5 X W^'^Mq), while any putative AGN in their centres 
(log(M,/M0) < 8.1) is probably smaller than those of lo- 
cal galaxies of similar mass (by a factor of more than 6) but 
with plenty of fuel to grow if it can tap the infered gas mass: 
the mass that a SMBH at our derived upper limit needs to 
accrete to reach the local SMBH-to-gas mass ratio is only 
about 1.3% of the total gas mass. 

The QSO RX J0941 shows strong X-ray emission (lu- 
minosity of 3 X 10** erg/s), high SFR (a few thousand so- 



lar masses per year) and strong ionized winds (Page et al 



2010). Within the evolutionary scheme outlined in the In- 
troduction, all these clues would mean that the SMBH is 
starting to push away the material responsible for the star 
formation and the SMBH growth. 

Finally, the properties of our SMGs are similar to those 



of similar objects in "blank field" observations (e.g. Borys et 



|al.|2005| [Alexander et al.|2008| whose properties have been 
obtained with similar methods) despite being in a signifi- 
cant overdensity. Similar properties of "spike" and "field" 
SMGs at z ~ 1.99 have also been reported by e.g. |Chap-| 



man et al. (20091. This indicates that the environment does 



not affect substantially the SMG properties. Galaxy inter- 
actions are commonly invoked to channel material to the 
internal regions of galaxies, triggering starbursts (Miho s fc] 



Hernquist 1994 1. We have mixed evidence: RX J0941 (850_1) 



and 450_4 appear isolated in the optical/NIR bands, while 
850_2 and 850_3 are in crowded immediate environments. 
Unfortunately, we lack high angular resolution data around 
our sources to be able to constraint tightly their interaction 
history. 



6 CONCLUSIONS 

Using X-ray to radio data, we have studied the immediate 
environment of RX J0941, a, z = 1.82 QSO with strong ion- 



ized winds visible in the UV and in X-rays ( Page et al. 2010 1, 



which shows as well both strong submm emission and an 
overdensity of submm galaxies (SMGs) around it |Stevens] 
et al.|2004 2010 1. Such high density regions are expected to 
form a galaxy cluster at ^ = (Kauffmann 1 996| . 

The 6 and 20cm radio data confirm the presence of 
all the submm sources and help to pinpoint their observed- 
frame optical-to-MIR counterparts. We have obtained pho- 
tometric estimates of their redshifts using hyperz to fit 
SWIRE galaxy templates I Rowan- Robinson et al.||2008 l to 
RiZJK and Spitzer-IRAC data. 

We found that four of the five unique submm sources are 
indeed associated with the QSO. The fifth source appears 
to be a background source, perhaps with an associated lobe 



robust against using different sets of templates ( [Bruzual fc 
Charlot|2003l|Rowan-Robinson et al.|20 08l) 



Under the assumption that a growing Super Massive 
Black Hole -SMBH- is hosted by our SMGs, we have used 
X-ray upper limits and the observed 24 pm fiuxes to estimate 
a very conservative upper limit to the X-ray luminosity (< 
10** erg/s) of such putative AGN (and hence a SMBH mass 
of log(M,/M©) < 8.1), even if Compton thick material is 
absorbing the direct AGN emission. 

We have used the rest-frame optical-to-radio starburst 



templates of Chary & Elbaz (2001 1 to estimate the star for 



mation rates (SFR) of the host galaxy of the QSO and the 
SMGs, obtaining individual rates of lOOOM0/yr or above 
(see Tablejl]). We have also estimated dust (~ 10® M©) and 
gas (~ 5 X IO^^Mq) masses, assuming parameters typical of 
the centres of SMGs. 

We have also estimated their stellar masses from their 
rest-frame A'-band luminosities, finding log(M,/M0) ~ 
11.5 ± 0.2, in line with other mass estimates for SMGs (e.g. 



Borys et al. 20051, but about 10 times larger than typical 



UV-selected star- forming galaxies at similar redshifts. 

An immediate conclusion (as also found for other sam- 
ples of SMGs) is that they are mature galaxies, with large 
stellar masses and an inferred reservoir of gas that would 
only allow a further 10% increase in mass at most. 

At the same time, their SMBH masses are a least factor 
of ~ 6 below that expected for their galaxy mass and the 



local SMBH-to-galaxy mass ratio ( Marconi & Hunt 2003 1 



This has also been found previously for other samples of 
SMGs, and is commonly interpreted as the growth of the 
SMBH lagging behind the galaxy growth, since the latter 
is already fully matured while the AGN still requires sub- 



stantial growth (Borys et al. 2005 Alexander et al. 2008 



Coppin et al. 20091. Accretion of a few percent of the in- 



fered gas mass on to the SMBH would be sufficient to reach 
the local mass ratio. 

Further observations are required to understand better 
the issues discussed here. Deeper X-ray data would allow to 
place tighter limits on, or to actually detect, AGN emission 
from the SMGs. Forthcoming Herschel SPIRE and PACS 
data will "fill in" the gaps in the FIR-submm SED, getting 
better estimates of the SFR and dust masses. High angular 
resolution millimetre observations would give more accurate 
estimates of the gas mass. Photometric redshifts and SED 
studies for the other overdensities of SMGs in our sample 
( Stevens et al.pOlO l will assess their relationship with the 
central QSOs and their evolutionary status, helping to un- 
derstand the role of these high density peaks in the forma- 
tion and evolution of galaxies. 
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APPENDIX A: APERTURE MATCHING OF 
THE PHOTOMETRIC POINTS 

Since each of our 5 optical-NIR images has a different seeing, 
we need to smooth all of them so that the effective seeing 
is the same as the one with the worst value {Z in our case). 
We have done this using GAUS in iraf. 

We have chosen a default aperture of 1.5" radius, which 
corresponds to ~98% enclosed flux for a 2D gaussian with 
FWHM equal to the effective seeing (or 96% for the em- 
pirical aperture correction explained below). There are 75 
sources closer than 3" to another source, so that the default 
aperture would assign counts from the same pixels to at least 
two sources. To avoid this, for each source we have chosen 
an aperture which is the minimum of 1.5" and half the dis- 
tance to the closest source (in steps of 0.1"). The minimum 
aperture found using this method is 0.6". 

The flux F of each source is determined (using phot 
within iraf) from the number of counts within the cho- 
sen aperture (after background subtraction), and the magni- 
tude from the usual — 2.51og_F. The default phot procedure 
to obtain the uncertainty in the magnitude involves deter- 
mining the mean value (msky) and the standard deviation 
(stdev) of the pixel values over the background area. This 
value can also be used to estimate the la limit for the flux 
of each source as ^msky x area (assuming Poisson statis- 
tics) where area is the area for the extraction of the source 
counts. 

Following this recipe produces both too low values for 
the magnitude errors and unrealistically low flux limits. We 
believe that this is because the reduction processes out- 
lined in Section [2] imply co-adding many individual images 
with re-scalings to match median or average pixel values 
over essentially arbitrary areas. Furthermore, since we have 
smoothed the images to have the same "effective" seeing, 
we have introduced correlations between the values of neigh- 
bouring pixels, and essentially averaging away the dispersion 
in the sky values. This is bound to alter the statistical char- 
acteristics of the counts in each flnal pixel, with respect to 
the naiVe Poisson assumption underlying most photometry 
packages (like phot). 

We have followed an alternative approach. For each im- 
age we have excluded an area of radius 2.5" around each 
source and, for each aperture radius, we have placed at ran- 
dom a large number (700) of circles of that size on the im- 
age, and we have measured the standard deviation (skydev) 
of the values of the fluxes within the simulated apertures, 
scaled to the total circle area. We have only kept simulated 
apertures in which more than 99% of their area was outside 
the excluded zones. The number of simulated apertures has 
been chosen not to oversample the un-excluded area. We 
have checked that reducing the fraction of good pixels to 
90%, increasing the exclusion radius to 3" or increasing the 
number of simulated apertures to 2000 does not change the 
obtained values. The error on the fluxes can then be esti- 
mated as -y/ F/epadu -I- skydev^ and the la upper limits as 
mo — 2.5 logj^Q(skydev), where epadu is the e~-per-ADU of 
the detector and mo are the zero point magnitudes given in 
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TablefT] The la limits obtained in this way are shown in Ta- 
ble |A|in terms of magnitudes. In order to get, for example, 
5(T limits, 2.51ogjQ(5) ~ 1.75 needs to be subtracted from 
those values. 

This problem is also present in the Spitzer images, so 
we have obtained skydev for them in the same way as for 
the optical-NIR images, using the source extraction radii 
given in Section |2.1[ Since those images were in units of 
MJy/sr, the errors on the fluxes obtained by SExtractor 
(assuming Poisson statistics) were underestimated, and very 
similar for all sources (since they were dominated by the 
variance of the pixel values in the background). We have 
used the FLUXCONV and GAIN keywords in the headers of the 
4.5 and 8 /im image^ and the exposure times to estimate 
the original number of counts in the image, and the corre- 
sponding errors. The total uncertainty in the image fluxes 
was obtained by adding to these errors in quadrature the 
total uncertainty in the determination of the values of those 
keywords (about 5%, see footnote) and the uncertainty in 
the sky values (skydev) . These fluxes and errors are the ones 
used for the empirical flux calibration described in |2.1[ After 
applying this calibration, the Icr sky standard deviations in 
the 4.5, 8 and 24 /im bands are 0.5, 2.1 and 52 /ijy, respec- 
tively. 

In order to obtain a good estimate of the stellar mass 
it is essential to have a good measurement of the rest-frame 
-ft'-band absolute magnitude, which would fall between the 
IRAC bands at the redshift of RX J0941. Hence it is essential 
to include those bands in the photometric redshift fits for the 
101 sources detected in either of those two MIR bands. But, 
since the Spitzer fluxes are corrected to "infinite" aperture, 
good estimates of this correction for the optical-NIR fluxes 
are needed as well, so that all the photometric points fed to 
the fit match. 

We have done this empirically, extracting for each 
source a radial profile in the smoothed i-band image between 
0.2" and half the distance to the nearest source (see above) 
and fitting to it both a single and a double gaussian profile, 
allowing thus for both "nuclear" and "galaxy" components. 
We have chosen for each source the best fit between those 
two using an F-test. The average value of the cr of a single 
gaussian profile (for the sources for which this was the best 
fit) is 0.59", slightly larger than the one used for matching 
the angular resolution of all optical-NIR images (0.54"), but 
essentially compatible with it. A correction factor was then 
calculated from the best fit profile for each source. 

The fiuxes/magnitudes of the 101 sources with 4.5 or 
8 ^m detections have been corrected with this factor be- 
fore the photometric redshift fitting, so the derived physical 
magnitudes are already for "infinite" aperture. For the re- 
maining 138 sources, this factor has been applied to these 
derived magnitudes after the photometric redshift fit. 

This paper has been typeset from a T^TjX/ M^TjX file prepared 
by the author. 



Table Al. la upper limits for each image and aperture (see 
App.|A]l 



Aperture 


R 


i 


Z 


J 


K 


(arcsec) 












0.6 


27.06 


27 Al 


26.26 


24.67 


23.44 


0.7 


26.76 


27.05 


26.02 


24.43 


23.08 


0.8 


26.55 


26.86 


25.85 


24.13 


22.86 


0.9 


26.28 


26.73 


25.64 


23.93 


22.65 


1.0 


26.15 


26.49 


25.49 


23.74 


22.51 


1.1 


25.95 


26.38 


25.40 


23.60 


22.34 


1.2 


25.80 


26.25 


25.21 


23.47 


22.21 


1.3 


25.68 


26.14 


25.10 


23.32 


22.08 


1.4 


25.51 


25.99 


25.01 


23.17 


21.99 


1.5 


25.40 


25.89 


24.89 


23.06 


21.91 
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